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We revisit the calculation of electroweak bremsstrahlung contributions to dark matter annihi- 
lation. Dark matter annihilation to leptons is necessarily accompanied by electroweak radiative 
corrections, in which a W or Z boson is also radiated. Significantly, while many dark matter models 
feature a helicity suppressed annihilation rate to fermions, bremsstrahlung process can remove this 
helicity suppression such that the branching ratios Br(£vW), Ht(£ + £~Z), and Yir(vvZ) dominate 
over ~Br(£ + £~) and Br(i'i'). We find this is most significant in the limit where the dark matter mass is 
nearly degenerate with the mass of the boson which mediates the annihilation process. Electroweak 
bremsstrahlung has important phenomenological consequences both for the magnitude of the total 
dark matter annihilation cross section and for the character of the astrophysical signals for indirect 
detection. Given that the W and Z gauge bosons decay dominantly via hadronic channels, it is 
impossible to produce final state leptons without accompanying protons, antiprotons, and gamma 
rays. 
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I. INTRODUCTION 

The importance of electroweak radiative corrections to 
dark matter annihilation has recently been recognized, 
and examined in a number of publications [lMldj. In 
a recent paper some of the present authors considered 
electroweak bremsstrahlung contributions to dark mat- 
ter annihilation, in models in which dark matter an- 
nihilation to a fermion-antifermion pair, XX ~ * f j \ is 
helicity suppressed There it was shown that W/Z 
bremsstrahlung lifts helicity suppressions, and can there- 
fore be the dominant DM annihilation mode. How- 
ever, some of the quantitative conclusions of [jj must 
be modified, as the explicit cross section calculation in 
therein was in error. The purpose of the present pa- 
per is to revisit and extend the calculation of the W/Z 
bremsstrahlung cross sections, and draw inferences from 
the result. The main inference is that the three body final 
state processes can still dominate the tree level process 
as claimed in [ij. We show herein that the claim finds 
support in the region where the parameter [i = m 2 /m 2 
is not too far from unity, with m r) and m x being the 
mass of the boson which mediates the annihilation pro- 
cess and the dark matter mass, respectively. This region 
of parameter space is reminiscent of the co-annihilation 
region in standard supersymmetric (SUSY) scenarios, al- 
though the present work can also be applied to models 
which are not in the SUSY framework. 

Let us parametrize the dark matter annihilation cross 
section in the usual way, 



a + bv , 



(1) 



where the constant a arises from s-wave annihilation 
while the constant b receives contributions from both s- 
and p-wave channels. Since the dark matter velocity in 



a galactic halo today is v ~ 10 _3 c, the p-wave term is 
strongly velocity suppressed. In order to have a large an- 
nihilation cross section in the Universe today, it is desir- 
able to have an unsuppressed a (s-wave) term. However, 
the s-wave annihilation of DM to a fermion-antifermion 
pair is helicity suppressed in a number of important and 
popular models. The most well known example is the 
annihilation of supersymmetric neutralinos to a fermion- 
antifermion pair. The circumstances under which helicity 
supressions do or do not arise were discussed in detail in 
Ref. Q. 

It has long been know that bremsstrahlung of photons 
can lift such a helicity suppression, leading to the re- 
sult that the cross section for XX ~^ f fl can dominate 
over that for XX ~> ff (llHl^ . However, the fact that 
radiation of a W or Z gauge boson would also lift a he- 
licity suppression had been overlooked until the work of 
Refs. In these scenarios for which the helicity sup- 

pression is removed, the dominant annihilation channels 
are the set of bremsstralung processes, namely 7, W and 
Z bremsstrahlung. (If the dark matter annihilates to col- 
ored fermions, radiation of gluons would also contribute). 
The phenomenology of W and Z bremsstrahlung is richer 
than that for photon bremsstrahlung alone. This is be- 
cause the W and Z bosons decay dominantly to hadronic 
final states, including antiprotons, for which interesting 
cosmic ray bounds exist. 



II. EXAMPLE OF SUPPRESSED 
ANNIHILATION 

To illustrate our arguments, we choose a simple ex- 
ample of the class of model under discussion. This is 
provided by the leptophilic model proposed in Ref. (l7l . 
Here the DM consists of a gauge-singlet Majorana 
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fermion x which annihilates to leptons via the SU(2)- 
invariant interaction term 



X + h.c. = f(v L T)° - i L V + )x + h.c 



(2) 

where / is a coupling constant, e is the 2x2 antisymmet- 
ric matrix, and (r] + , -q ) is a new SU{2) doublet scalar. 
In this model, DM annihilation to fermions is mediated 
by t and u channel exchange of the rj fields. 

An identical coupling occurs in supersymmetry if we 
identify x with a neutralino and r\ with a sfermion dou- 
blet. In fact, the implementation of supersymmetric 
photinos as dark matter by H. Goldberg provided the 
first explicit calculation of s-wave suppressed Majorana 
dark matter annihilation to a fermion pair [l9j . There- 
fore, much of what we discuss below is also relevant for 
neutralino annihilation to fermions via the exchange of 
sfermions. However, the class of models for which the 
2 — > 2 annnihilation is helicity suppressed is more gen- 
eral than the class of supersymmetric models. 

The cross section for the 2 — ► 2 process XX 
e + e~ or vv is given by 



/V i 

247rm2 (1 



M) 4 



(3) 



where m; ~ and m„± = m n o have been assumed, and 
fi = m 2 /m 2 . The suppressions discussed above are ap- 
parent in Eq. Q. The helicity suppressed s-wave term 
is absent in the mi =0 limit, and thus only the v - 
suppressed term remains. 



III. LIFTING THE SUPPRESSION WITH 
ELECTROWEAK BREMSSTRAHLUNG 

A. W-strahlung cross section 

We shall take the limit mi ~ and assume that m j; ± = 
m„o. The matrix elements for the six diagrams of Fig. [1] 
are given by 



M a = i 



■f 2 9 1 1 



V2 qf h 



2 -f„ TTl 2 



Mh = i 



M c = -i 



x(v(k 2 )P L v{p 2 ))(u( Pl )^P L ^u(h))e^, (4) 

f 2 g i i 

y/2 qf ui - m 2 
x(«(fc 1 )P LW (p 2 ))(S(p 1 ) 7 ^P L ^ u (fc 2 ))e« (5) 

f 2 9 1 1 



\/2 ql t 2 - ™ 2 



M d = -i 



x(v(k 2 )P L ^v(p 2 ))(u(p 1 )P R u(k 1 ))e^, (6) 
f 2 9 1 1 



V2 ql u 2 



x{{v{k 1 )P L ^v{p2)){u{jp l )P R u{k 2 ))e^, (7) 



Me = -i 



f 2 9 



M f = -i 



V2h-m 2 t' 3 - m 2 

x{{v{k 2 )P L v{p 2 )){u{ Pl )P R u{ki)) 
x((h-Pi) + (ki-pi-Q)reQ, (8) 

f 2 9 1 1 



m 2 u' 3 



y/2u 3 ..„„ „ 3 ... v 

x{{v{k l )P L v(p2)){u(p l )P R u{k 2 )) 

x((fc 2 -pi) + (fc 2 - P i-Q))^, (9) 

where we define the usual helicity projectors Pr/l = 
i ( 1 ± 75), and the Mandelstam variables 

ti = (fci - qi) 2 , 
t 2 = (fci - pi) 2 = <3, 
ui = (fc 2 - <ji) 2 , 
u 2 = (fc 2 - pi) 2 



«3, 

(k 1 - Pl -Q) 2 , 
u' 3 = {k 1 -p 2 f = (k 2 - Pl -Q) 2 . 



t' 3 = (k 2 - P2 ) 2 



The vertex factors used in the matrix elements are as 
follows: the IvW vertex has an -^= r y fJ - P L e^ , the xvl ver ~ 
tex is i/Pt, and the coupling between the W~~ and the 
rj + —jfl is taken to be of the form —ig(p + p')/y2 from 



Ref . [20| • Fierz transformed versions of these matrix ele- 
ments, and some insight gained from them, are collected 
in Appendix [A] 

We have explicitly checked the gauge invariance of our 
set of Feynman diagrams. Writing the matrix element as 

M =M»e®, 



the Ward identity 



Q M M M = 0, 



(10) 



(11) 



is satisfied for the sum of the diagrams. The Ward iden- 
tity takes the same form as for photon bremsstrahlung 
provided we take the lepton masses to be zero, since 
the axial vector current is conserved in this limit. Note 
that diagrams (a) + (c) + (e) form a gauge invariant sub- 
set, as do (b)+(d) + (f). The full amplitude is the sum 
of the partial amplitudes, properly weighted by a minus 
sign when two fermions are interchanged. Thus we have 
M = {M a + M c + M e ) -(M b + M d + M /). 

In performing the sum over spins and polarizations, we 
note the standard polarization sum, 



pol. 



'w 



(12) 



can be replaced with —g^u alone. The Ward identity 
of Eq. lfTT]) ensures the second term in Eq. (|12|) does not 
contribute once the contributions from all diagrams are 
summed (and squared). 

In addition, we find that the longitudinal polarization 
of the W also does not contribute to the s-wave ampli- 
tude, i.e. 



(13) 
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(a) 



«-(9l) 



"(pi) X(fcl) 




"(Pi) X(*l) 



(<0 




ffe) 



"(Pi) X(k\) 



W-(Q) 



e + (P2> x(*2) 



W-(Q) 



(e) 



(f) 




"(pi! 



e + ( P2 ) 
"(Pi) 



W-(Q) 



e + (p 2 ) 



FIG. 1. The i-channel ((a),(c), and (e)) and it-channel ((b), (d) and (f)) Feynman diagrams for \X ~ ► e, + vW~ . Note that t- 
and M-channel amplitudes are simply related by the ki <-¥ k^ interchange symmetry. All fermion momenta in the diagrams flow 
with the arrow except P2 and ga, with q\ = p\ + Q, qi = f>2 + Q- 



The W boson behaves as a massive transverse photon, 
with just two transverse polarizations contributing. As a 
consequence, our calculation of W bremsstrahlung must 
reduce to the known results for photon bremsstrahlung 
in the mw limit, modulo coupling constants. Below 
we will show that this happens. 

The thermally-averaged cross section is given by 



v da = — / — 
2s I 4 



\M\ 2 dLips 3 



(14) 



spin, pol. 



where the \ arises from averaging over the spins of the 



initial \ pair, v = y 1 is the mean dark matter 

velocity, as well as the dark matter single-particle velocity 
in the center of mass fram^E and dLips n represents n- 
body Lorentz invariant phase space. 

We calculate the cross section for W emission follow- 
ing the procedure outlined above, with the integration 
over phase space performed according to the method de- 
scribed in Ref. [1]. We expand in powers of the DM 
velocity, v, keeping only the leading order (v ) contribu- 
tion. As expected, we have an unsuppressed cross section 
given by 



awf 4 
2567r 2 m 2 , 




-2Li 



-2 In 



2m2( M -l) 



-2Li 2 



In 



2(74(^ + 1)2-771^) 



4(m 2 x (ji + l) 2 -m 2 w ) 



(4/i + 3) m 2 w (4m 2 .(/j + l)(4/x + 3) - (m 2 w - 4m 2 x )(fi - 3)) 



(M+l) 



16m^(/i + l) 



m 2 w (4m^(/i + l) 4 - 2m 2 w m 2 x (fx + l)(p + 3) - m^Li - 1)) 



In 



2m%(n- 1) 



l) 3 (m 2 (/i 



In 



'w 



4m 2 



(p-l)(2m 2 x (fi + l)-m 2 w ) 



In 1 - 



m 2 (p + l) 5 



2m 2 x (fj, + 1) - m 2 w J Am\(p + l) 3 (4m 2 /z - m 2 w ) (m 2 x (p + l) 2 - m 2 w ) 



x (Am x {p + 1) 4 (4/j + 1) - m x m 2 w (p + l) 2 (3/i(/i + 6) + 7) + 2m£m&, + 4) + 1) - m 6 w ) \ (15) 



1 Informative discussions of the meaning of v are given in |2 ill , and 



the inclusion of thermal averaging is covered in [22 
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where aw = 9 /(47r) . The Spence function (or "dilogarithm") is defined as Li2(z) = — j Q ^ In |1 — C| = J2k 

If we take the limit mw —> and replace aw with 2a em , then Eq. (|15p reproduces the cross section for 
bremsstrahlung of photons, namelj0 



ij 



128tt 2 to2 



(/i + 1) 



2- - In 2 
6 



1 



2/y 



-2Li 2 



1 



2/t 



4/i + 3 4^ 2 - 3a* - 1 . /> - 1 



1 



2// 



In 



/' 



1 



(16) 




FIG. 2. The ratio R = v a( X X 
as a function of fi = (m,,/m x ) 2 , 



->■ e + i/W~ 
for m x = 



)/u<t{xX 
300 GeV. 



->• e ' e J 
We have 



used v = 10 c, appropriate for the Galactic halo. 




nur [GeV] 



FIG. 3. The ratio R = va( X X 
as a function of the DM mass m 
used v = 10 _3 c, appropriate for the Galactic halo. 



e vW )/v cr(xx e e ) 
x , for fj, = 1.2 GeV. We have 



The successful recovery of the photon bremsstrahlung 
result in the massless W limit provides a checlo 
on the rather complicated expression for massive W 



2 Note that Eq.2. of Ref. [3 is larger by an overall factor of two, 
and also has the opposite sign for the (l+fi)[...] term, while Eq. 1. 
of Ref. GJ is consistent with our results. 

3 A related work 23] appeared on the arXiv nearly simultaneously 
with ours. In this related work there appears analytic expres- 



> 
b 




m x [GeV] 

FIG. 4. The cross sections for XX ~* e + vW~ (red) and XX ~^ 
e + e~7 (blue), for fi = 1.2 and coupling / = 1. For large DM 
mass, the cross sections differ by a factor of 1/(2 sin 2 8\y) = 
2.17 while for m x comparable to mw the W bremsstrahlung 
cross section is suppressed by phase space effects. 



bremsstrahlung given above in Eq. (|15|1. 

Since we are working in the limits v — and rrif = 0, 
the nonzero results in Eqs. (fT5]) and ([TBI imply that 
the leading terms are neither helicity nor velocity sup- 
pressed. Not clear from the mathematical expressions 
is the sensible fact that the cross sections fall monoton- 
ically with increasing (or fx). This monotonic fall 
is shown in Fig. [31 where we plot the ratio of the W- 
strahlung cross section to that of the lowest order process, 



R = va(xx -> e + vW )/va(xX 



-). The lowest 



order process itself falls as /j,~ 2 , so the VF-strahlung pro- 
cess falls as \i~ A . This latter dependence is expected 
for processes with two propagators each off-shell by 1/fi, 
thereby signaling leading order cancellations among Fig. 
1 diagrams (a)-(d). 

Importantly, the effectiveness of the W^strahlung pro- 
cesses in lifting suppression of the annihilation rate is 
evident in Fig. [2] The ratio is maximized for /j, close to 
1, where m x and are nearly degenerate. However, the 
TT-strahlung process dominates over the tree level anni- 
hilation even if a mild hierarchy between m x and m v is 
assumed. The ratio exceeds 100 for \i < 2. 

Fig. [3] illustrates that the ratio R is insensitive to the 
DM mass, except for low m x where the W mass signifi- 



sions for the Mz , Mw = limits of the cross-section which we 
calculate, thereby providing another calculational check. 
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cantly impacts phase space. From the figure one gleans 
that for m x > 3mw, the ratio R is already near to its 
asymptotic value. Incidentally, the asymptotic value may 
be obtained analytically by dividing Eq. Q16p with Eq. ^ 
and rescaling a em with 0^/2. 

In Fig. U we compare the IF-strahlung cross section 
with that for photon bremsstrahlung. For high dark mat- 
ter masses where the W mass is negligible, the two cross 
sections are identical except for the overall normalization, 
which is higher by factor of 1/(2 sin 2 Oy/) = 2.17 for W- 
strahlung. For lower DM mass, the available phase space 
is reduced due to W mass effects, thus the VF-strahlung 
cross section falls below that for photons. This can be 
seen in Fig. 0] for m x < 150 GeV (this number is fairly 
insensitive to fj,) . Another factor of two is gained for W- 
strahlung when the W + mode is added to the W~ mode 
shown here. 

Nominally, the correct dark matter energy fraction is 
obtained for early-Universe thermal decoupling with an 
annihilation cross section of 3 x 10~ 26 cm 3 /s. It is seen 



in Fig. d] that the IF-strahlung mode falls 2-3 orders of 
magnitude below this value. Note that at the time of 
dark matter freeze-out in the early Universe, the velocity 
suppression of the p-wave contribution is not as severe as 
it is for late-Universe annihilation. Hence, radiative W- 
strahlung with its natural suppression factor aw/^ is 
probably not the dominant annihilation mode responsible 
for early-Universe decoupling of Majorana dark matter. 

B. W and Lepton Spectra 

To obtain the energy spectrum of the W, we compute 
the differential cross section in terms of Ew by making 
the transformation 

^o S (O q )^ {s _-^ m2w) dE w . (17) 

The energy spectrum of the the primary leptons is cal- 
culated in similar fashion. We find 



da 

dx 



IV 



awf A 
128tt 2 to2 



(1 - x w ) + 




(p + l)(p + l-2x w ) + 



(p + i- x w y 



((// + 1)(M + 1 - 2x w ) + Zjjf) ^ /(1 + 1-XW + ^x 2 w - m 2 w /ml 



2(p + 1 - xw) z 



fj,+ 1 — xw — \ x 



- rnwl m2 :, 



(18) 



da otwf A 1 
'dx~i ~ 5127r 2 m 2 (n-l + 2xif 



4(l-a; I ) 2 -4a;i(A» + l) + 3(/i + l) 2 



x In 



2m 2 (^ + l)(l-a; ; )-™ 2 



xi (4m 2 (1 - xi) - mlv) 



(2m 2 x (n + 1 - 2xi) - m 2 w ) (1 - x t ) I (2m|(l - xi)(p + 1) - m 2 w ) (1 - Xl f 



(1 - Xl ) 2 (4(1 - Xl f - Xl (p + 1) + 3(a* + l) 2 ) + - x t ) { Xi (jjl + 11) - 4(/i + 3)) ~xi^ 



.(19) 



The W spectrum per XX ~~ ^ svW event is given in 
Fig. [5J We use the scaling variable xw = Eyy /m x , and 
plot dN/dxw = ( ) a *+"W- . The kinematic range 

of x w is (1 + with the lower limit corre- 

x 

spending to a IF produced at rest, and the upper limit 
corresponding to parallel lepton momenta balancing the 
opposite W momentum. As evident in Fig. [SJ the W bo- 
son spectrum has a broad energy distribution, including 
a significant high energy component. 

For the lepton energy spectrum, shown in Fig. |6l the 

2 

range of the scaling variable xi = Ei/m x is [0, 1 — j^r]. 

Both limits arise when one lepton has zero energy and the 
other is produced back-to-back with the W. Note that 
this spectrum is valid for either the e + or the v from the 



annihilation XX ~~ ^ & + vW , and for either e or f from 
the annihilation XX ~~ ^ &~vW + . 



C. Z Emission 

Consider the process producing the vvZ final state. 
The cross sections for the Z-strahlung processes are re- 
lated to those for TF-strahlung in a simple way: The am- 
plitudes producing vvZ arise from the same six graphs 
of Fig[TJ where e, W and rj + are replaced everywhere 
by v and Z and 770, respectively. The calculation of the 
amplitudes, and their interferences, proceeds in an iden- 
tical fashion. After making the replacement raw fnz, 
the cross section for the annihilation process XX vvZ 
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0.4 0.6 0.8 

FIG. 5. The W spectrum per \X ~~ ^ evW annihilation for 
m x = 300 GeV and // = 1.2. 




0.4 0.6 

xi=Ei/m x 

FIG. 6. The primary lepton spectrum per XX e.vW anni- 
hilation, for m x = 300 GeV and fi = 1.2. 



differs from that for XX ~ * e + vW by only an overall 
normalization factor, 



1 



VU V 9Z 



(2 cos 2 e w ) 

0.65 x va e + uW 



x V (J e +uW- 



(20) 



Consider now the e + e~ Z final state. Again, the ampli- 
tudes arise from the same six basic graphs of FigfT] Since 
only the left-handed leptons couple to the dark matter via 
the SU(2) doublet rj, only the left handed component of 
e~ participates in the interaction with the Z . Therefore, 
the couplings of the charged leptons to Z and W take 
the same form, up to a normalization constant. We thus 
find 



VCJe+e-Z 



2 (sin 2 6 w -\)' 

cos 2 9w 
0.19 x va e + uW - 



x va e + U Y/- 



mw—rTnz 



(21) 



Adding the four contributions to VF/Zstrahlung, we find 

VOWyZ-strahlung = 2.84 X Va e + vW - . (22) 

IV. DISCUSSION AND CONCLUSIONS 

There are clear advantages and disadvantages of seek- 
ing photon- versus VF/Z-bremsstrahlung as an indirect 
signature of dark matter. With photon bremsstrahlung, 
the photon itself is easily detected. It's energy spec- 
trum may then be readily compared to model predic- 
tions. With lU-strahlung, it is the decay products of the 
W decay which must be sought. Their spectra are less 
attributable to the model of dark matter annihilation. 
However, the total rate for IU/Z-strahlung exceeds that 
of photon-strahlung. Photons couple with strength e, 
W's couple with strength = e/(V2sin9w), and 

Z's couple to neutrinos with strength g / (2 cos 6w) — 
e/(2 cosOw sin#w)- Therefore in the high energy limit 
where the W and Z masses can be neglected, we expect 



1 



Ve+vW 



2 sin 2 9 W 



&e+e-~f = 2.17CT e + e - 7 . 



(23) 



So, in the high energy limit where m x > 300GeV >> 
mw, the total cross section becomes 



fbrem, total — &e+vW- + <J ue-\V+ 
+OvvZ + &e+e-Z + 
= 7.16<7 e + e - 7 . 



(24) 



Furthermore, the varied decay products of the W/Z al- 
low more multi-messenger experiments to engage in the 
dark matter search. Charged leptons, protons and an- 
tiprotons, neutrinos, and even deuterons are expected, 
at calculable rates and with predictable spectra. Impor- 
tantly, hadronic decay products are unavoidable, despite 
a purely leptonic tree-level annihilation. The tens of mil- 
lions of Z events produced at CERN's e + e _ collider show 
in detail what the branching fractions and spectra are for 
each kind of decay product. In a forthcoming article [9] 
we explore the favorable prospects for using VF-strahlung 
decay products as indirect signatures for dark matter. 

The lifting of the helicity suppression is most signifi- 
cant in the limit where the mass of the boson mediat- 
ing dark matter annihilation does not greatly exceed the 
mass of the dark matter particle. This is true both for 
photon bremsstrahling and for PF/Z-bremsstrahlung. In 
this limit, we find the three body final state annihilation 
channels can significantly dominate over two body anni- 
hilation channels. The region of parameter space where x 
and r\ are approximately degenerate is of great interest in 
many models, since it coincides with the co-annihilation 
region where both XX an d XV annihilations are impor- 
tant in determining the relic dark matter density at the 
time of freezeout in the early Universe, often a favored 
parameter region in SUSY scenarios. 
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Appendix A: Fierz transformed matrix elements 

Upon Fierz transforming; (for standard 2 — ¥ 2 Fierz 
identities, see e.g., (2~3 . [25j|) the matrix elements of 
Eqn.g}-© we find 



[Pl] (fPLfJ = ~ Tr[P L T c i Q P Lil Y B ] [Y B ] [T c ) 



M a = 



M b = 



*9f 2 



1 1 



-(v(k 2 )-f a PRu(k 1 ) 



y/2qf h-m 2 

^/2 q f ui - m 2 2 V 
x(u{Pi)r^l a PLv{p 2 ))tf, 



(Al) 



(A2) 



-wf 2 1 1 



-(v(k 2 )j a PR.u(ki) 



y/2q% ti--"' 1 

—in f 2 1 1 / 
M d = —^r T - {v{k 2 )~f a P L u(ki) 



(A3) 



Mi 



s/2q% u 2 -m 2 2 
x(u(Pih a ^PLv(p 2 ))e^ 
-igf 2 (2fci - 2 Pl - QY 



2^2 (t' 3 - m 2 )(t 3 - m 2 ) 
x(u(p 1 ) 1 a P L v( P 2))e^ 
-igf 2 {2k 2 -2 Pl -QY 

2V2 K - m 2 )(u 3 - m 2 ) 
x(u( Pl )~f a P L v(p 2 ))eQ. 



(A4) 
v{k 2 )-f a P R u{ki)) 

(A5) 

{v{k 2 )~i a P L u(ki)) 
(A6) 



Alternatively, we may apply a chiral version of the 
Fierz transform (discussed in detail in Ref. to trans- 
form Eqs.Q-©. After a bit of algebra we get a pleasant 
factorized form for the bilinear currents. We show de- 
tails for the first one, and then summarize the results for 
current products of the other matrix elements. 

The current product in amplitude M a of Eq. Q is 

(v(k 2 )P L v(p 2 )) (fi(pi)^P L ^u(A;i)) . (A7) 
We write this current product in Takahashi notation [25] 



- - A Tr[P L j a f Q P L ^ } (P R >f] [P Lla ) . 

(A8) 

In going from the first equality to the second, we insert 
the only values for T c and Tb allowed by the helicity 
projectors in the string of gamma matrices. Finally, we 
may invert the sequence in the trace, and remove the 
Takahashi notation to write the result as 



1 



A Tr[P R ^P L f lin a ] x 

(«(pi)iW«(pa)) (v(k 2 )P Lla u(ki)) ■ (A9) 

Amplitude M& is computed in a similar way. In ad- 
dition, it is useful to use the identity for a Majorana 
current 

(v(ki)P L ~i a u(k 2 )) = {v(k 2 )P R ~t a u(ki)) [Majorana] 

(A10) 

to put the final result in a form similar to that for am- 
plitude M a . The other current products are reduced in 
a similar fashion. The final result for the product of cur- 
rents after Fierzing is 

i(z;(fc 2 )P 7 Q U (fc 1 )) (u( Pl )P R ^v( P2 )) (All) 
f Tr[P R f$ llpla ] , for M a , M b 
x I Tr[P L fi^ a } , for M c , M d 
{ 2 g a /3 , 



for M e , M f 



In addition, the unspecified projector P in the first com- 
mon factor is Pj, for amplitudes M a , M c , M e , and P R 
for the amplitudes Mb, Md, Mf derived from the crossed 
graphs. 

What can we learn from this exercise? For graphs M e 
and Mf the Fierzed currents are the same as in the 2 — > 2 
case. This fact is not surprising since in these graphs the 
internal W emission does not perturb the form of the 
currents and their product. However, for the other four 
graphs with W emission occurring on a fermion leg, the 
form of the current product is quite different from the 
2^2 case. With 2 — !• 3 scattering, the Lorentz in- 
dex of each current need not contract directly with the 
other. Referring to Table 1 of Ref. [1], one sees that 
unsuppressed Majorana annihilation amplitudes become 
possible for the axial vector combination (757 ] [757V )j 
and for the vector combination (7 s ] [7V ) , providing the 
trace post-factors in Eq. (| Al 1|) do not vanish. These 
combinations are at the heart of the unsuppression which 
we have presented in this paper. (The role of amplitudes 
M e and Mf is to cancel gauge non-invariant contributions 
from graphs M a -Md-) See also Ref. (2(| for a comprehen- 
sive discussion of enhanced/suppressed DM annihilation 
modes. 

Also, for m 2 >> t, u, the non-current factors in am- 
plitudes M a and M& are the same, as are the non-current 
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factors in amplitudes M c and M<j. Then the subtraction 
of one from the other leads to a pure axial vector coupling 
in the Majorana current. This in term leads to an effec- 
tively pure axial vector coupling in the final state lepton 
current. This effective axial vector-axial vector coupling 



of currents was advertised earlier. However, for values 
of t and u which are non-negligible when compared to 
mi, there is some residual vector coupling. In this more 
complicated case, it is probably best to directly calculate 
rates without Fierzing the currents. Such is the course 
followed in the main text of this paper. 
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